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Abstract: The Mandelbrot set is an extremely well-known mathematical 
object that can be described in a quite simple way but has very interesting and 

04 , non-trivial properties. This paper surveys some results that are known concern- 

ing the (non-)computability of the set. It considers two models of decidability 

r \ , over the reals (which are treated much more thoroughly and technically in [I], 

j2], and 0] among others), two over the computable reals (the Russian school 
and hypercomputation) and a model over the rationals. 
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1 Introduction 



In theoretical computer science it is no surprise that the halting problem for 
Turing machines is the favourite target for solution by non-conventional models 
of computation. However, the decidability of sets of reals and the computability 
f^ \ of functions in ordinary real analysis is a topic of great interest to the broader 

mathematical community and a potential area of application for so-called hyper- 
computation. This paper surveys some of the most important results and gives 
an extremely simple example of the application of accelerated Turing machines 
to the question of decidability of the Mandelbrot set. This very amusing prob- 
lem was raised by Roger Penrose jS] and — like many good questions — implies 
considerable work on the definitions, in this case what exactly is meant by de- 
5_^j ■ cidable for a subset of the plane. 

In 1979 Benoit Mandelbrot used a computer to plot 1 a beautiful approxima- 
tion of the subset 

M = {c e C | for all n > 1, |/"(0)| < 2} where f c (x) = x 2 + c (1) 

of the complex plane C (where /" denotes n-th iteration of /) . This set was orig- 
inally described by Pierre Fatou in 1905 but after the appearance of a colourful 
plot of the set in Mandelbrot's book and — especially — on the cover of Sci- 
entific American and in an accompanying column in August 1985, the Man- 
delbrot set has become one of the greatest celebrities of popular mathematics. 

'Department of Decision Sciences, University of South Africa (Pretoria), PO Box 392, 
UNISA, 0003, Republic of South Africa, potgiph@unisa.ac.za. 
1 Mandelbrot actually plotted a mirror image of M. 




Figure 1: Membership candidates for the Mandelbrot set 



Plots can easily be generated using a plethora of software freely available, in- 
cluding normal IATj^X code 2 . The relatively concise but non-optimised Octave 3 
code snippet 



# For an nxn grid 

# Number of iterations 

# Set up grid 

# Initial value on grid 

# Iterate the mapping 

# Plot monochrome, absolute 

# value of 2.1 is escape 

suffices, for example, to plot the Mandelbrot image in Figure ^ 

The image in Figure^shows an expanse of background (white) points which 
were shown during the execution of the code to have left the closed disk of 



n=1000; 

m=50; 

c=meshgrid(linspace(-2,2,n) )\ 

+i*meshgrid(linspace(2,-2,n) ) ' 
x=zeros(n,n) ; 
for i=l:m 

x=x. ~2+c; 
endf or 
imagesc(min(abs(x) ,2.1)) 



^ http://www.thole.org/manfred/apf el/apf el.tex [accessed 2005-12-30] by Manfred 
Thole. The Mandelbrot set is approximated by the white area at the centre of this plot 
generated by a correctly compiled apfel.tex 

3 Octave is a free and open-source high-level language for numerical computation with 
implementations on many platforms and largely compatible with MATLAB®. See 
http : //www . octave . org/ 




Figure 2: The Mandelbrot plot as we too often see it 



radius 2 — therefore not to belong to the Mandelbrot set — and grey and black 
points which are candidates for membership of M. In this plot the points c have 
a lighter colour when they have approached relatively closely to the boundary 
of the disk after exactly 50 iterations of the map / c (0). This shading scheme 
emphasises that the usual routines for generating Mandelbrot plots, like this 
one, clearly can identify only elements of the complement of M and not of the 
set itself. Scientists are well-aware this fact (that the complement of M is only 
known to be recursively enumerable) but it is usually not overemphasised in 
the more popular writing. Incidentally, M is compact and simply connected 8 
which one cannot, and should not expect to, see in Figure ^ Figure 2 shows a 
plot of the Mandelbrot set candidates as one more usually sees it. 

The central part is, of course, the same as the set of grey or black points in 
Figure n — those points who have not (yet) escaped the disk after m iterations 
in the Octave code fragment 

n=1000; # For an nxn grid 

m=60; # Number of iterations 

c=meshgrid(linspace(-2,2,n) )\ # Set up grid 

+i*meshgrid(linspace(2,-2,n) ) ' ; 
x=zeros(n,n) ;y=ones(n,n) ; # Initial values on grid 

# y counts number of iterations 

# without escape from the disk 
for i=l:m 

x=x.~2+c: # Usual iteration 



y=y .+( .5+sign(3-abs(x) ) . /2) ;# Add one if still in radius 3 
x=x. * (min(abs(x) ,3)\ # Scale back points far away to 
. /(abs(x)+ ! abs(x) ) ) ; # speed up and avoid overflow 

endf or 

imshow(imagesc(y) ,rand(m+l ,3) ) # Plot y with random colours 

which now keeps track of when a point first leaves the disk (if at all) using 
the matrix y. The union of the outer coloured bands represents R \ M on the 
computational grid, or what we could discover of it in the number of iterations 
executed. Although pretty, these bands show the stages of construction of M's 
complement and disguise the fact that we know about M not much more than 
that it is somewhere in the central coloured section! Incidentally, John Ewing 
and Glenn Schober have made [5] laborious numerical estimates of the area 4 of 
the Mandelbrot set using two different plotting methods and arrived at answers 
of respectively 1.52 and 1.72. This is for a set which everybody believes that 
they have seen! 

In a very readable article ^J| Lenore Blum has described the gulf between 
numerical analysis and computer science. Blum calls the classical theory of 
(Turing) computability 

fundamentally inadequate for providing such a foundation for 
modern scientific computation, in which most algorithms with ori- 
gins in Newton, Euler, Gauss, et al.-are real number algorithms 

as a partial justification of the Blum-Shub-Smale model of computation over 
the real numbers (mentioned in 13. II below). The present author would be in- 
clined to the opinion that 'modern scientific computation' is rather inadequately 
founded in the classical theory of computability and that the logical and founda- 
tional problems which arise in this regard (some of which are illustrated in this 
paper) could be much deeper and connected to our model of the real numbers 
and their representation; and to the finite/infinite duality which is drilled into 
mathematical recruits in a dingy room in the Hilbert Grand Hotel. 

2 M in the Recursive Realm 

First, let us consider only computable points in the Mandelbrot set M. By this 
we mean points that can be finitely specified and communicated in a consistent 
way, which one can take to mean represented by a program for computing approx- 
imating rationals. This is an intuitive idea used in all versions of computable 
analysis and which can be readily grasped. Fix a universal Turing machine U 
as well as a (recursive, surjective) encoding <f> : N —* {0, l}xfjxN with 

cj):n^ {(j}i(ri),<j>2{ri),(j)z(ri)) 

so that every rational number is of the form 



4 A tighter estimate has recently been computed by Yuval Fisher and Jay Hill |10l . 



for some n. If x € R then we say that a; belongs to the set R c of computable 
reals if there exists a program 5 i x for [/ so that, if /, x is the function computed 
by i x then 
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for every m G N. Penrose's question in this section should be: 

Given a program-description of a computable complex number, can 
we algorithmically determine whether it belongs to M, or not? 

In other words, does there exist a (partial) function G :C N x N — > N, com- 
putable in some sense, such that whenever i x and i y are programs for x, y G M c 
respectively then 

G(i x ,i y ) = %M(x + iy). (2) 

This simply means that G can be used to determine membership of M for the 
recursive points in the plane. Identify C with R 2 in the usual way and consider 
M c = MnR 2 ., the computable points in the plane that belong to the Mandelbrot 
set. The unsuspecting reader of popular scientific literature could reasonably 
assume that given a full description of a point (x,y), by a program pair (i x , i y ) 
one supposes there exists a procedure for deciding membership of M c . The 
existence of such a G has been implicitly suggested to the general public for two 
decades by the pretty pictures we have had our computers draw, but it turns 
out that we do not have such a procedure in classical computability theory at 
all. 

2.1 Markov computability 

In the Russian school of constructive mathematics, pioneered by Andrei Markov, 
a function / : M c — > K c is considered to be computable (or, constructive) if 
there exists a Turing machine computable G : N — » N such that whenever i x 
is a [/-program for x then G(i x ) is a [/-program for f(x) — and f(x) is defined 
if and only if G(i x ) is. A set is computable or decidable in this setting when 
its characteristic function is computable. However, every computable function 
is continuous J2| an d therefore the only computable sets will be closed (and 
open). This is a rather ironic development as, for example, the computable 
interval [0, 1] fl M c will not be a computable set in this sense, being closed 
but not open in K c . The same is true of the unit circle and the disk in the 
plane and in both cases the difficulty springs from the fact that there is not 
general procedure that, given i x and i y , can decide whether x — y or not. The 
Mandelbrot set is, however, closed in R (and hence closed in R c ), so could it 
be computable in this sense? No, since —2 G M c and —2 is a cluster point in 
M c of M c \ M c and hence M c is not open in R c . Consequently its characteristic 
function cannot be Markov computable 6 . 

5 We shall assume that the program and input for U (actually both simply inputs) are 
natural numbers. 

6 Since the unit circle is also not computable here, this should not come as a great surprise. 



Although apparently intuitive, the notion of computable set used in this 
subsection is clearly very, very bad from the point of view of real analysis. Nev- 
ertheless it corresponds in some sense exactly to what a programmer would 
regard as computable: given a procedure (subroutine) for an x, having a pro- 
gram that outputs 1 in finite time if x S M c and zero otherwise. The situation 
here is in contrast to that of subsets of the natural numbers, for which the no- 
tion of decidability is very natural and well-established (pace the entire field of 
'super- Turing' hypercomputation) . 

2.2 Zeno Machine Computability 

Consider what may be called a Zeno machine (ZM) |KS| or Accelerated Turing 
Machine ^J^l- With this kind of speed-up of the computing device, one can 
solve the halting problem for Turing machines in finite time. Many hypercom- 
putational schemes tend to be proposals for somehow accomplishing infinitely 
many computational steps in a finite time (see also |16j . for example). Without 
loss of generality, we shall assume a ZM to be identical to a Turing machine with 
one input tape, one output tape and a storage tape except that the ZM takes \ 
hour to execute the first transition, ^ hour for the second, i hour for the third 
etc. After one hour the ZM will have finished its operation and one will perhaps 
find the answer to some tantalising question on the output tape. On a putative 
ZM one could implement an Octave interpreter that would execute the code 

n=1000; # For an nxn grid 

c=meshgrid(linspace(-2,2,n) )\ # Set up grid 

+i*meshgrid(linspace(2,-2,n) ) ' ; 

x=zeros(n,n) ; # Initial value on grid 

do 

x=x.~2+c; # Usual iteration 

x=x. * (min(abs (x) ,3)\ # Scale back points far away to 

. /(abs(x)+ ! abs(x) ) ) ; # speed up and avoid overflow 

# and infinite values 

until (1==0) # Repeat a lot 

imagesc(min(abs(x) ,2 . 1) ) # Plot x, 2.1 counts as escape 

in finite time. This would provide an exact plot of the Mandelbrot set on 
the grid points! With a small modification (adding a procedure/subroutine for 
approximating the computable real) a ZM can decide membership of M c for 
any computable real. For example, the Octave code immediately above can be 
rewritten for an (ordinary, Turing computable) Octave function c(i) (instead of 
a matrix c) where c(i) gives a rational approximation of c to within 2 -4 and 
a scalar x, initially zero. In the i-th iteration of the loop we then recompute 7 
x using c(i). A similar calculation could be done relative to an oracle for the 
halting problem. The ZM as used here does not apparently present any of 

7 This would require recomputing the previous iterates, of course, but there are only finitely 
many to do each time. 



the problems with respect to defining the terminal configuration of the devices 
described in J3| since the matrix x is always bounded. This stability is however 
illusory: in executing the code described here one needs to continually reset 
x back to zero and therefore the variable x will, for every x € R c \ M c , have 
values alternatively and with absolute value 2.1 arbitrarily close to the end of 
execution time 8 . 

2.3 A Rational Refuge? 

Is there some relief from these problems if we restrict our attention to the 
points with rational coordinates only? Consider again the (recursive, surjective) 
encoding <fi : N — > {0, 1} x N 2 of the rational numbers used earlier. A set A C Q 
can be defined as computable whenever a Turing machine computable function 
F : N — ► N exists such that 

^l^-^A) =XAO(j). 

In this sense, now, the rational points on the unit circle do constitute a com- 
putable set since the condition x 2 + y 2 = 1 can be checked by a Turing machine 
for rational x and y. Is M n Q computable in this sense? 

This kind of computability over the rationals is very different from Markov- 
computability over the computable reals. Consider for example the function 
/ : Q —> {0, 1} such that f(q) = 1 if the reduced improper fraction representa- 
tion of q has an even denominator and f(q) = otherwise. This function is not 
continuous on Q and therefore not the restriction of a Markov-computable func- 
tion to the rationals. It is therefore conceivable 9 that a Turing machine could 
compute the characteristic function of M (~l Q with respect to the representation 
4> of the rationals, -F|</,-i(MnQ)- 

The rational points are perhaps a bad basis for developing a general theory 
of computability of subsets of R n however. It could say nothing much about 
the curve x 3 + y 3 = 1. In fact, decidability with respect to the rationals suffers 
from a general failure to take the boundary into account (as in the example 
below). Nevertheless, for connected and compact sets with non-empty interior, 
computability in this sense seems a relatively natural and quite desirable prop- 
erty. It would, for instance, allow one to plot the set with a computer using 
test points on a rational grid. Consider also that the set {(x,y) 6Q 2 | y > e x } 
is computable with respect to the representation (f>: if ei(q,m) is a computable 
function approximation of e q from below and e2(q,m) from above such that 
lim„ woo ei(q, m) = e q then an enumeration of the values ei(g, m) and e2(q, rn) 
for m = 1, 2, . . . will after finitely many steps reveal whether any given rational 
lies below e q or above it (since e q is irrational for all rational q ^ and the case 
q = can be checked separately). 

8 Thomson's Lamp showing the way... 

9 Although it strikes the present author as unlikely that M will be computable in this sense, 
a proof is called for. 



3 M in Real Space 

Let us return to the standard real numbers and consider M as a subset of the 
standard plane. In pursuing an answer to Penrose's question, the formulation 
of an appropriate notion of decidability of subsets is again required. Most of 
the work in this regard has its roots in the Polish school of computable analysis, 
starting with Andrzej Gregorczyk and Daniel Lacombe. The Blum-Shub-Smale 
system on the other hand, has a rather algebraic flavour. 

3.1 The Blum-Smale Result 

Blum, Mike Shub and Steven Smale (BSS) have introduced [3 ^] a model of 
computation over arbitrary commutative rings which is based on machines that 
operate using the elements of the ring R in lieu of a finite alphabet. There exist 
universal machines in this model. In the case where the ring is Z2 the classical 
computability theory is recovered. BSS-computable functions over R are func- 
tions computed by such a machine and they call a set computable or decidable 
whenever its characteristic function is BSS-computable. Blum and Smale have 
shown that the Mandelbrot set is NOT computable in this framework 17 which 
is at least a partial answer to Penrose's question. However, Vasco Brattka has 
shown |3] that in the BSS scheme over the field of standard real numbers the 
set {(x,y) G R 2 | y > e x } is not computable either. Brattka's result reflects 
unfavourably on the claim that BSS computability provides a natural notion of 
decidable set — at least, for real analysis. 

3.2 Computable Real Analysis 

Computable real analysis in the Polish school, as developed in ^B] and else- 
where, is based on the definition of a function as computable whenever it maps 
every computable sequence of points (in R") to a computable sequence of points 
and has a recursive modulus of continuity (defined on {1, |, |, j, . . .}) on every 
compact subset. The appropriate definition [3j in this context is for a set A 
to be computable whenever its distance function cLa '■ R" — ► R is a computable 
function in this sense. Computable (in the classical sense) subsets of N, viewed 
as subsets of R, remain computable in this sense pQ and the notion is there- 
fore a true generalisation of the notion of classical Turing computability. Peter 
Hcrtling has recently shown [I] that if the Mandelbrot set M is locally con- 
nected 10 , then its distance function has to be computable and hence Penrose's 
question would be answered in the affirmative. 

As described in [3|, in this notion of computability the closed decidable 
subsets of R™ are exactly those sets which can in principle be plotted with 

10 Actually, Hertling proved a stronger result: that the hyperbolic conjecture (which would 
be implied by local connectedness) would be sufficient to prove imply computability of the 
distance function. 



arbitrary accuracy on a computer screen 11 . It is therefore at this point not yet 
known whether the Mandelbrot set can be accurately drawn by a computer! 

4 Conclusion 

The table below summarises the results and the two open problems mentioned 
in this survey. 

Circle y > e x M'brot 
Markov-computability over M c x x x 

Blum-Shub-Smale over R / x x 

Turing-computability over Q ■/ ■/ ? 

Computable analysis (Brattka e.a.) / ■/ ? 

Zeno-computability over M c •/•/•/ 

Although the classical computability of the rational Mandelbrot set M (~l Q 2 
is an obvious question which one would cautiously expect to be answered in 
the negative, the author is not aware of a current result implying this. Among 
the models of decidability of sets in R" , the approach 13.21 studied by Brattka, 
Weihrauch e.a. appears the most reasonable and a demonstration of the com- 
putability of M in this setting would be extremely interesting both in itself 
as well providing strong support for the intuitiveness of their approach. In the 
model over Q the Mandelbrot set is of course decidable with respect to an oracle 
for the halting problem but it seems unlikely the converse is true, so if M n Q 2 
is not Turing-decidable then it could be interesting 12 to study problems that 
can be solved relative to an oracle for MdQ 2 (or an M-oracle in any model 
in which M is not decidable) . Until such time as at least one of the question 
marks in the table have been decided, a Zeno machine (or any hypercomputing 
model capable of solving the halting problem for Turing machines) remains — 
alas! — the best way of imagining that we can actually decide membership of the 
Mandelbrot set. 
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